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Abstract 

The baby Skyrme model is a (2+l)-dimensional analogue of the Skyrme model, 
in which baryons are described by topological solitons. In this paper we introduce a 
version of the baby Skyrme model in which the global 0(3) symmetry is broken to 
the dihedral group Djy. It is found that the single soliton in this theory is composed 
of N partons, that are topologically confined. The case N = 3 is studied in some 
detail and multi-soliton solutions are computed. These soliton solutions are related 
to polyiamonds, which are plane figures composed of equilateral triangles joined by 
common edges. It is shown that the solitons may be viewed as pieces of a doubly 
periodic soliton lattice. An alternative model with D3 symmetry is also introduced, 
which has an exact explicit soliton lattice solution. Soliton solutions are computed and 
compared in the two D3 theories. Some comments are made regarding the extension 
of these ideas to the Skyrme model. 
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1 Introduction 



The Skyrme model [1] is a nonlinear field theory in which baryons are described by topolog- 
ical solitons, called Skyrmions. The model has been obtained from quantum chromodynamics 
(QCD) as a low energy effective field theory in the limit in which the number of colours, N, 
is large [2]. More recently, the Skyrme model has been derived from string theory, in the 
context of holographic QCD, again in the large N limit [3]. 

The number of colours, N, appears in the Skyrme model only as a coefficient of the 
Wess-Zumino term. This plays an important role in the quantization of Skyrmions, but at 
the classical level the soliton solutions are blind to the value of N; as the Wess-Zumino term 
does not contribute to the classical energy. An interesting issue is whether it is possible to 
incorporate an effective small value of iV at the level of the classical soliton solution. 

The baby Skyrme model [4] is a (2+l)-dimensional analogue of the Skyrme model, that 
has proved to be a useful testing ground for the study of Skyrmions. Soliton solutions of the 
baby Skyrme and related models are also of interest in their own right, within the context 
of condensed matter physics [5, 6], where direct experimental observations can be made. 

In this paper we introduce a version of the baby Skyrme model in which the global 0(3) 
symmetry is broken to the dihedral group D^. It is shown that this reproduces some key 
features expected of Skyrmions in a toy model associated with a small number of colours. 
In particular, it is found that the single soliton in this theory is composed of N partons, 
that are topologically confined. The case N = 3 is studied in some detail and multi-soliton 
solutions are computed. The model admits a variety of stable multi-solitons that take the 
form of polyiamonds; which are plane figures composed of equilateral triangles joined by 
common edges. Doubly periodic soliton lattices are also computed and it is shown that the 
polyiamond solitons may be viewed as pieces of a soliton lattice. 

An alternative model with D$ symmetry is also introduced and studied. This model has 
the property that an exact soliton lattice solution is given explicitly in terms of a Weierstrass 
elliptic function. Soliton solutions are computed and compared in the two D 3 theories. 

2 The broken baby Skyrme model 

The field of the baby Skyrme model is a three-component unit vector <j> = (<f>i, </> 2 , ^3). 
In this paper we are concerned with static solitons, hence the theory may be defined by its 
static energy, which takes the form 

E = J (^drf ■ d t <t> + x dj4>) ■ (drf x dj4>) + V^j d 2 x, (2.1) 

where V(<f>) is a potential and k is a constant. In the standard baby Skyrme model [4] the 
potential is taken to be 

\/ = m 2 (l-0 3 ), (2.2) 

which is the analogue of the conventional pion mass term in the Skyrme model [7]. The 
constant m gives the mass of the fields 4>i an d 02, associated with elementary excitations 
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around the unique vacuum <fi = (0, 0, 1). 

Finite energy requires that the field takes the vacuum value at all points at spatial infinity, 
<fi(oc) = (0,0, 1). This compactification means that topologically </> is a map between two- 
spheres, with an associated integer winding number B e Z = 7r 2 (S' 2 ). This topological charge 
(or soliton number) is the analogue of the baryon number in the Skyrme model and may be 
calculated as 

B = -— J <t> ■ (&0 x d 2 (f)) d 2 x. (2.3) 

An application of Derrick's theorem [8] reveals that the scale of a soliton in the baby Skyrme 
model is determined by the ratio ^Jn/m. 

The first two terms in the energy (2.1) are invariant under the global 0(3) symmetry, 
4> i — y 0<p for O G 0(3), but this is broken by the potential (2.2) to an 0(2) symmetry acting 
on the first two components 0i, 2 . The single soliton takes advantage of this symmetry and 
is axially symmetric [4]. 

Other choices for the potential have been investigated [9, 10], but in these examples there 
is an unbroken 0(2) symmetry and the B = 1 soliton is axially symmetric. A novel situation 
was considered recently [11] using the easy plane potential V = ^m 2 ^. Again this leaves an 
unbroken 0(2) symmetry, but in this case the choice of a vacuum value at spatial infinity, 
for example 4>{oo) = (0,0, 1), distinguishes a point on the orbit of the unbroken symmetry 
and further breaks the symmetry to the dihedral group D 2 . In this case the single soliton 
is not axially symmetric, but turns out to be composed of two constituents. The work on 
easy plane baby Skyrmions provided some inspiration for the model proposed in the current 
paper, but it should be stressed that the two theories are quite different, as described shortly. 

The only previous work we are aware of in which the potential has only a discrete sym- 
metry is the work of Ward [12], in which the symmetry is broken to the dihedral group D 2 
by the choice V = \m 2 (l — — 4>\). In this theory the single soliton is also composed of 
two constituents. We shall discuss this theory in more detail in section 5, where we introduce 
a generalization to a model with D 3 symmetry and compare the results with those obtained 
for the theory of main concern in this paper. 

The theory introduced in the present paper has the potential 

\/ = m 2 |l-(0 1 + # 2 ) iV | 2 (l-03), (2.4) 

where N > 2 is an integer parameter of the model. We shall refer to the theory with this 
potential as the broken baby Skyrme model with N colours, to use a suggestive notation. 
The global 0(3) symmetry is broken by this potential to the dihedral group Dn, generated 
by the rotation (0i + ifa) h-> (0i + i0 2 )e l27r/,JV and the reflection (0i, 2 , 03) > (0i 5 ~ 02, 03)- 
In the context of symmetry groups in three dimensions, this pyramidal symmetry group 
is often denoted Cjy v , but as an abstract group it is isomorphic to D N , which is a more 
convenient notation for the later application to planar symmetries. 

The potential (2.4) has N + 1 vacua on the two-sphere. The vacuum at the north pole 
<fi = (0,0,1), will be the chosen vacuum at spatial infinity, with the remaining iV vacua 
lying on the equatorial circle 3 = at the iVth roots of unity. Note that this choice of 



3 



boundary condition does not break any further symmetries, in contrast to the situation for 
the easy plane potential. Another important difference between the two theories concerns 
the masses of the elementary excitations of the fields 0i and 02- To quadratic order in <pi and 
02, the broken potential (2.4) agrees with the standard baby Skyrme potential (2.2), hence 
both fields 0i and 2 have mass m. However, for the easy plane potential clearly only the 
0i field has mass m and the 2 field is massless. Hopefully, this brief discussion has served 
to highlight some of the important differences between the two models, making it clear that 
different phenomena should be expected, even in the case N = 2 where both theories have 
the same unbroken symmetry group D 2 . 

Ward's potential [12] with D 2 symmetry has more in common with the two colour broken 
baby Skyrme model, as both fields 0i and 02 have the same mass. However, there is an 
additional inversion symmetry <fi i— >■ —(f), with an associated extra vacuum <p = (0,0,-1), 
which leads to some qualitative differences in the soliton solutions. 

3 Solitons and polyiamonds 

Whether in two or three dimensions, a symmetry of a soliton in a Skyrme model refers to 
an equivariance in which a spatial rotation or reflection can be compensated by the action 
of the global symmetry of the theory. Clearly the maximal symmetry possible in the broken 
baby Skyrme model with N colours is dihedral symmetry D^. Even the B = 1 single soliton 
cannot be axially symmetric, though the expectation is that it has the maximal symmetry 
Dm. Numerical solutions presented in this section confirm this expectation and reveal that 
the single soliton is composed of N constituents, which we refer to as partons, given that in 
the Skyrme model the single soliton describes the proton. 

Each parton carries baryon number B = 1/N, associated with a winding that covers 
this fraction of the target two-sphere. Partons are topologically confined, since finite energy 
requires that the total baryon number is integer-valued, hence an equal number of partons 
of each colour. Here the term colour may be used to enumerate the N segments of the 
target two-sphere obtained by drawing great semi-circles from the north pole to the south 
pole that pass through each of the N equatorial vacua. Of course a parton and an anti- 
parton (associated with a winding of the same segment with the opposite orientation) is an 
allowed combination, corresponding to elementary excitations of the 0i and 02 fields with 
zero baryon number, being the baby Skyrme analogue of pions. 

From now on we shall concentrate on the most physically relevant case of N = 3 colours, 
regarding the theory as a lower-dimensional analogue of the Skyrme model. The generic 
values k — m — 1 are taken for the parameters of the theory. 

To numerically construct soliton solutions an energy minimizing gradient flow algorithm is 
applied to the energy (2.1). Spatial derivatives are approximated using fourth-order accurate 
finite difference approximations with a lattice spacing Ax = 0.05 and 25 1 2 grid points. At 
the boundary of the grid the field is fixed to the vacuum value <f> = (0, 0, 1). 

A field with topological charge B is given by 



4> = (sin / cos(B9) , sin / sm(B9) , cos /) , 
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where r and 9 are polar coordinates in the plane with f(r) any monotonically decreasing 
radial profile function such that /(0) = 7r and / vanishes at the boundary of the grid. 

The field (3.1) has dihedral symmetry D 3 b, with the spatial rotation 9 i— > 9 + 2n/(3B) 
being compensated by the global transformation (0i + 202) i— >■ (0i + i02)e _l27r//3 , and the 
spatial reflection 9^—9 compensated by the global reflection (01,02,03) >— > (01, —02,03)- 
Note that either a spatial reflection or a global reflection alone changes the sign of the 
topological charge, but the combination of the two leaves it unchanged. 

Using the field (3.1) as an initial condition in the numerical energy minimization code 
yields a charge B solution with D 3 b symmetry. The energy density of the B = 1 soliton is 
displayed as a contour plot in Figure 1.1, which clearly shows the three constituent partons. 
A plot of the topological charge density (the integrand in (2.3)) has a similar structure. The 
energy of this D3 symmetric B = 1 soliton is E = 34.79, and is listed in Table 1. For all the 
numerical results presented in this paper, the topological charge computed using the lattice 
version of (2.3) is integer-valued to five significant figures, which provides an indication of 
the accuracy expected in the numerical computations. 




Figure 1: Energy density contour plots for solitons with charge B. The top row displays 
stable solitons for B = 1, 2, 3 and the bottom row shows three different stable solitons with 
B = 4. 



The D 3 b symmetric solution is composed of 3B partons located at the vertices of a 
regular 35-gon. The solutions with B = 2 and B = 3 are displayed in Figures 2.2a and 
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2a 




Figure 2: Energy density contour plots for solitons with charge B. 
2a) B = 2 hexagon; 3a) B = 3 nonagon; 3b) B = 3 triquetra. 
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Table 1: The energy per soliton E/B and symmetry group G of solitons with topological 
charge B < 4. 

2.3a respectively. For B > 2 such solutions are unstable to perturbations that break the 
dihedral symmetry. The B = 2 hexagonal solution is stable with an energy per soliton of 
E/B = 33.04. A sufficiently large symmetry breaking perturbation can convert this solution 
into the additional stable D 2 symmetric soliton displayed in Figure 1.2. The energy per 
soliton of this solution is E/B = 33.06, so the two different B = 2 solutions have energies 
that differ by an amount comparable to our expected numerical accuracy, and we are unable 
to make a confident statement about which has the lower energy. It is clear from Figure 1.2 
that this B = 2 soliton is constructed from two B = 1 solitons, with a relative spatial 
rotation of 180°. 

As mentioned earlier, the asymptotic fields of a soliton in the broken baby Skyrme model 
have the same form as those in the standard baby Skyrme model; hence the results on 
asymptotic interactions derived in [4] can be transfered to the current theory (see Appendix 
A for a discussion of asymptotic forces for a general potential). In particular, the leading 
order result shows that two single solitons are maximally attractive if one is rotated relative 
to the other through an angle of 180°. As the single soliton in the broken baby Skyrme 
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Figure 3: The six possible polyiamonds constructed from at most four triangles. 



model is not axially symmetric, then beyond leading order there will be a contribution that 
differentiates between the relative orientation of the two triangles of partons. The result 
presented in Figure 1.2 displays the optimal orientation between the two triangles, which 
can be confirmed by using an initial condition consisting of two well-separated single solitons 
with a generic initial orientation. 

It is instructive to represent the single soliton, as in Figure 3.1, as a triangle with three 
coloured dots to denote the three peaks in the energy density associated with the three 
segments of the target two-sphere. With this representation the B = 2 soliton in Figure 1.2 
corresponds to the diamond in Figure 3.2, in which the two triangles share a common edge 
with adjacent dots being different colours. 

This triangle representation suggests that multi-solitons will correspond to polyiamonds. 
A polyiamond [13] is a plane figure composed of identical equilateral triangles joined by 
common edges, so that no two triangles overlap. For three triangles there is a unique 
triamond, shown in Figure 3.3. Given a polyiamond, and a colour assignment to any one 
of the triangles, there is a unique colouring of the triangles satisfying the rule of different 
adjacent colours. 

The initial condition (3.1) with B = 3 yields the unstable D$ symmetric nonagon dis- 
played in Figure 2.3a with an energy per soliton of E/B = 33.68. A symmetry breaking 
perturbation leads to the B = 3 soliton displayed in Figure 1.3 with E/B = 32.83. This 
confirms the predicted triamond form of Figure 3.3, which has only a reflection symmetry, 
D\, and no rotational symmetry. This B = 3 soliton can also be obtained from an initial 
condition consisting of three well-separated single solitons with appropriate positions and 
orientations. A stable B = 3 local energy minimum has also been computed that is not 
of the polyiamond form. It has an energy per soliton of E/B = 33.00 and is presented in 
Figure 2.3b. This triquetra solution has D% symmetry and is formed from three triangles 
that share vertices but not edges. 
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With four triangles there are the three tetriamonds shown in Figures 3.4a, 3.4b and 
3.4c. Suitable initial conditions, using four well-separated single solitons with appropriate 
positions and orientations, leads to B — 4 solitons associated with each of these tetriamonds. 
Energy density plots for all three solitons are displayed in the bottom row of Figure 1 and 
the associated energies are listed in Table 1. All three solutions appear to be stable and their 
energies are very close to each other. This can be understood from the fact that all pairs of 
triangles that share a common edge are in a maximally attractive orientation. Furthermore, 
to transform from one configuration to another requires the breaking of an attractive bond 
followed by a change of orientation and a repositioning of a triangle. 

The number of polyiamonds grows rapidly with the number of triangles, and the expecta- 
tion is that there will be a soliton associated with each of these. There are five pentiamonds 
(including the first asymmetric configuration with neither a rotation nor a reflection sym- 
metry) and twelve hexiamonds. Already for B = 9 there are 160 polyiamonds, so it is likely 
to be a computationally intensive task to numerically compute all the solitons for values of 
B larger than those considered in this paper. However, it could be a worthwhile exercise 
that might lead to an interesting energy function on the space of polyiamonds. Based on 
the result for B — 4, it could be that the linear arrangement of triangles is minimal for all 
values of B. This would have some similarities with the standard baby Skyrme model, where 
soliton chains are found to be the minimal energy configurations [14]. 

Note that the B = 3 triquetra solution in Figure 2.3b is related to the B = 4 tetriamond 
solution in Figure 1.4b by removing the central triangle. A local energy minimum with 
B = 4 has been computed that is not of the polyiamond form but is related to the triquetra 
solution by the addition of a triangle (in the polyiamond manner) on an outside edge of the 
triquetra, rather than in the middle. It is expected that a variety of similar local energy 
minima exist for all B > 4, based on extending the triquetra solution in this manner. 



The polyiamond solutions in Figure 1 suggest the existence of a doubly periodic triangular 
lattice, associated with a tiling of the plane by equalateral triangles. To study a lattice with 
triangular symmetry it is sufficient to restrict to the fundamental torus T 2 with a 60° angle 
in the (x±,X2) plane. It is useful to note that several identities may be derived for such a 
lattice to be a critical point of the energy (2.1) (see Appendix B). 

The first identity is a standard virial relation that follows from an application of Derrick's 
theorem [8] and is a requirement of criticality under a rescaling of the lattice 



The two remaining identities follow from variations of the lattice associated with a stretch of 
one of the fundamental periods and a variation of the angle of the fundamental torus from 



4 Soliton lattices 



k 



f 




(4.2) 



60° 




(4.3) 
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and 

J r d l( j>-d 2 cj>d 2 x = 0. (4.4) 

For computational purposes it is convenient to work with a rectangular torus of the form 
(xi,X2) G [0, x [0, a/3L], containing two copies of the fundamental torus. Numerical 
simulations can then be performed in the rectangular torus with periodic boundary conditions 
in both the x\ and x 2 directions. An initial value is chosen for the length L of the torus and 
the energy minimized using the numerical methods described in the previous section. The 
virial relation (4.2) is then evaluated and used to predict an improved estimate for the torus 
length L. This procedure is iterated to convergence and finally the two remaining identities 
(4.3) and (4.4) are checked. For the results presented in this paper the rectangular lattice 
contains 52 x 90 grid points and all three identities are satisfied to an accuracy of better 
than 1%. 

The fundamental torus of a triangular tiling contains two triangles, hence the rectangular 
torus contains a field with topological charge B = 4. An initial condition is provided by 
setting B = 4 in the ansatz (3.1) and using a radial profile function with a compact support 
inside the rectangular torus. A periodic perturbation is applied to this initial condition, to 
break any reflection symmetries. The result of the numerical minimization is displayed in 
the left-hand side of Figure 4. The minimizing torus length is L = 5.38, with an energy per 
soliton for this lattice equal to E/B = 31.83. This value is consistent as a limit of the finite 
B polyiamond energies presented in Table 1. 




Figure 4: Energy density contour plots for two different soliton lattices. On the left is a 
triangular lattice of single solitons and on the right is a hexagonal lattice of double solitons. 
The single soliton lattice has a slightly lower energy per soliton. 

There is a second soliton lattice, based on the hexagonal B = 2 soliton of Figure 2.2a. 
This lattice is obtained if the initial condition is not subjected to a perturbation that breaks 
the left-right reflection symmetry. This lattice is displayed in the right-hand side of Figure 4. 
In this case the minimizing torus length is slightly reduced at L = 5.10 and the energy per 
soliton is a little greater at E/B = 32.43. The fact that the energy of the double soliton 
lattice is larger than that of the single soliton lattice agrees with the preferred polyiamond 
form of the minimal energy solitons with B > 2. 
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All the polyiamond solutions in Figure 1 may be viewed as finite B pieces cut from 
the soliton lattice presented in the left-hand side of Figure 4. This explains why different 
solutions with the same value of B have very similar energies. The difference is an edge 
effect associated with the way in which a finite B piece is cut from the infinite lattice. This 
supports the expectation that stable soliton solutions exist for each polyiamond and suggests 
that the energy differences between solutions will decrease as B increases. 



5 An alternative D 3 model 

As mentioned earlier, in [12] Ward introduced a potential with D 2 symmetry. It has the 
novel feature that an exact explicit solution exists for the soliton lattice, with the single 
soliton composed of two partons. The lattice has square symmetry and consists of half 
solitons. In this section, we apply Ward's analysis to construct a potential with _D 3 symmetry, 
for which there is an exact explicit triangular lattice solution. We then investigate the solitons 
of this alternative D$ model, pointing out some of the differences and similarities with our 
earlier three colour theory. 

Consider the energy (2.1) defined on a torus T 2 , with an arbitrary potential V. Let b 
denote the topological charge density, the integrand in (2.3), 

6 = — f (^x9 2 # (5.1) 

47T 

The inequality 

(Annb - VW) 2 > (5.2) 
leads to the following lower bound for the energy on the torus 

E > Air^B + k J^ 2 bV2Vd 2 x^j, (5.3) 

where, without loss of generality, we have restricted to the situation where B > 0. 
Following [12], the inequality (5.3) can be written as a Bogomolny bound 

E/B > 4tt + Kp, (5.4) 

where (5 is the constant 

J2 , 



(3= / V2Vd 2 S, (5.5) 
J s 

with d 2 S the standard area element on the target two-sphere, normalized to 4ir. 

To investigate the bound further, it is useful to introduce complex coordinates on both 
the domain and target by defining z = x\ + ix 2 and W — (0i + «0 2 )/(l — 03). A necessary 
condition to attain the bound (5.4) is that W is a meromorphic function of z. With this 
assumption the topological charge density (5.1) becomes 



b = 



1 


dW 2 


tt(1 + \W\ 2 ) 2 


dz 



(5.6) 
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The bound is then attained if the inequality (5.2) becomes an equality, which requires 



V 



8k 2 


dW 


(1 + 


W 


2)4 


dz 



(5.7) 



This defines a suitable potential if dW/dz is an explicit function of W. As the domain is 
a torus, the most elegant choice is to take W to be proportional to a Weierstrass elliptic 
function. For simplicity, we set W(z) = p(z), where the Weierstrass function is defined in 
terms of the invariants g 2 ,g3 by 

(%) 2 = ~ 92P ~ to- (5-8) 

Substituting this choice into (5.7) yields the family of potentials 

T/ 8k 2 \AW 3 - g 2 W - g,\ 2 



Ward studied the case g% = 0, when the potential (5.9) has D 2 symmetry (W h- >■ iW), with 
an exact lattice solution of half solitons and square symmetry [12]. In this paper we consider 
the theory with g 2 = 0, so that there is D 3 symmetry (W >->■ e l27T ^ 3 W). Different values 
of gz are related by a scaling symmetry, which we use to set gs = 4. The potential of our 
alternative D 3 model is therefore given by 

128k 2 I1 — W 3 \ 2 

V = (1 + 1^2)4 = 16ft2 ( 1 - + 30^ + Witt - 0?), (5.10) 

where in the final expression the potential is written in terms of the field </>. 

As in the previous sections, from now on we set k — 1. The vacua of the potential (5.10) 
are the same as those of the three colour broken baby Skyrme model (2.4). Considering 
elementary excitations around the vacuum = (0,0,1), reveals that the fields <f>i and 02 
both have mass m — 8. This is not an important difference, and is simply a consequence of 
our choice of scaling for the exact elliptic function solution, but does mean that the scale of 
the solitons will be smaller in the alternative D 3 model and the energies higher. 

A significant difference between the two D 3 theories concerns the relative values of the 
potential at particular points of interest on the target two-sphere. These four points are 
<fi = (0, 0, —1), defining the centre of the soliton, and the three points on the equator midway 
between the equatorial vacua, such as </> = (— 1,0,0). In the broken baby Skyrme model the 
potential at a midway point is greater than the potential at the soliton centre, producing a 
split of the single soliton energy density into three peaks. However, for the alternative D 3 
theory this situation is reversed, so the energy density is not split but is simply stretched 
into a triangular deformation. In this sense the alternative model does not reveal the parton 
constituents in the same way as the broken baby Skyrme model. However, as we now 
describe, there are a number of similarities between the solitons of the two theories. 

The energy density of the exact elliptic function solution is displayed in Figure 5. To aid 
comparison with the previous lattice solutions, the energy density is plotted for a rectangular 
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Figure 5: An energy density contour plot of the lattice of single solitons given by an exact 
elliptic function solution of the alternative D 3 model. 



torus containing two copies of the fundamental torus. The length of the rectangular torus is 
equal to the real period of the elliptic function, hence L = u\ = r(|)r(|)/(2v / 37r) = 2.428... 

A numerical evaluation of the integral (5.5) yields /3 = 89.30 and hence from (5.4) an 
energy per soliton of E/B = 101.86. Note the similarity between Figure 5 and the left-hand 
side of Figure 4. Both are triangular lattices of single solitons, with the main difference being 
that the energy density of the single soliton is not split into three peaks in the alternative 
D 3 theory, as discussed above. The similarity between the minimal energy lattices in the 
two theories suggests that solitons in the alternative D 3 theory may also be related to 
polyiamonds. To investigate this issue, the same numerical code used in section 3 is applied 
to the problem, with the only change being a reduction in the lattice spacing to Ax = 0.02, 
reflecting the smaller scale of the solitons. 
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Table 2: The energy per soliton E/B and symmetry group G of solitons with topological 
charge B < 4 in the alternative D 3 model. 

Energy density contour plots are displayed in Figure 6 for some solitons with B < 4. The 
similarity to Figure 1 for the broken baby Skyrme model is obvious, as is the polyiamond form 
of the solutions. The symmetry and energy per soliton for each of these solutions is presented 
in Table 2 and the results are consistent with the limiting lattice value of E/B = 101.86. 

One difference between the two theories concerns the non-polyiamond solutions, such as 
the B = 2 hexagon solution shown in Figure 2.2a. In the alternative D 3 model a B = 2 
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Figure 6: Energy density contour plots for solitons of charge B in the alternative D% model. 
The top row displays stable solitons for B = 1, 2, 3 and the bottom row shows three different 
stable solitons with B = 4. 



hexagon appears to be unstable, with an energy per soliton of E/B = 106.33, which is larger 
than that of the diamond solution E/B = 105.80. 

6 Conclusion 

A version of the baby Skyrme model has been introduced in which the global 0(3) sym- 
metry is broken by the potential to a dihedral symmetry Dm, with the result that the single 
soliton is composed of N topologically confined partons. Multi-solitons have been computed 
in the N = 3 theory and shown to be related to polyiamonds, with a form consistent with 
cutting pieces from a doubly periodic soliton lattice. It would be worth extending the re- 
sults in this paper to larger soliton numbers, to verify that the polyiamonds correspondence 
continues. It might also be amusing to investigate the N = 4 theory and determine whether 
the multi-solitons are related to polyominoes. 

The obvious extension of the work in this paper is to the (3+l)-dimensional Skyrme 
model. It is straightforward to construct an analogous symmetry breaking potential in the 
Skyrme model, though the physical consequences of isospin symmetry breaking must be 
carefully considered. The potential term in the Skyrme model does not play as crucial a 



13 



role as in the baby Skyrme model, so the influence of a symmetry breaking potential is not 
expected to be as dramatic. However, even the inclusion of the traditional pion mass term 
does produce qualitative differences in the structure of multi-Skyrmions, for sufficiently large 
baryon numbers [15, 16], so some new features should appear. 



Appendix A: long range inter-lump forces 

In this appendix we present an analysis of the long-range forces between solitons in the 
baby Skyrme model with a general smooth potential V : S 2 — > [0, oo). Let v e l /_1 (0) be 
the vacuum value of <fi (so lim^oo cp = v). By definition, this is a (possibly degenerate) 
minimum of V, so the Hessian of V at v is well-defined, and has non-negative eigenvalues 
> 0, which we choose to order so that fa\ < /j, 2 ,. (Recall that the Hessian of V at v is 
the unique self-adjoint linear map J v : T V S 2 — > T V S 2 such that 



d 2 V(d>(t)) 



dt 2 



= X-J V X, (6.1) 

=o 



where <f>(t) is any curve in S 2 with 4>(0) = v, 4>(0) = X.) The case fa\ = was treated in 
[11], so we shall concentrate on the case where fa\ > 0. Let e±, e 2 be the corresponding unit 
eigenvectors, oriented so that e 2 = v x e±. At large r one expects the fields of a soliton to be 
well approximated by a solution of the linearization of the field equation about v, namely 

<f>(r, 9) — v - m^Ki (far) cos Qe x - -77-^1 (far) cos 6e 2 H (6.2) 

where K v denotes the modified Bessel function of the second kind and (71,(72 are unknown 
real constants which will receive a physical interpretation shortly. We define linearized fields 
Xii X2 so that <fi = v + xi^i + X2 e 2 and observe [4] that (6.2) corresponds to the solution of 
the linearized model, with Lagrangian density 

C = E (^Xi^Xi - \Ax\ + (6-3) 

in the presence of external point sources 

«i = qid x S(-x), k 2 = q 2 d y 5(yL). (6.4) 

Asymptotically, the soliton coincides with the fields induced by scalar dipoles of moment 
qi pointing along the x-axis and q 2 along the y-axis, inducing fields of mass fa and fa 
respectively. This is the soliton in standard orientation and position. If the soliton is rotated 
through an angle a and translated to a G M 2 , it corresponds to the composite point source 

«i = Qi • V5(x - a), k 2 = q 2 • V5(x - a), (6.5) 

where qi = q\ (cos a, sin a) and q 2 = g 2 (— sina, cosa). The interaction energy experienced 
by two solitons placed at a and a with orientations a, a is expected to coincide asymptotically 
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as R = |a — a| — >■ oo with that of their corresponding point sources Kj, /?j interacting via the 
linear theory (6.3), 

U = - J (kiXi + K2X2) d 2 x, (6.6) 

where x% denotes the field induced by A lengthy but straightforward calculation yields 
the formula 

U = T,%^ { K °(* R )«M + \ii - otti] ) , (6-7) 

where || and _L represent the components of q relative to the orthonormal basis n = (a— a)/R 
and rr 1 . 

The form of U depends strongly on whether fii and /i 2 are equal. If fii < /i 2 , the 
expression (6.7) is dominated at large R by its first term 

U=£K ( t i 1 R)qli + ... (6.8) 

which predicts that the force between two solitons is maximally attractive when the dominant 
dipoles qi, qi are anti-aligned along the line joining the soliton centres (qi = — qi = ±q± n). 
The neglected terms in (6.7) may become significant at intermediate range, particularly if 
H 2 2 — [if is small. 

If Hi = H2 and qi = q 2 , then whatever the orientations of the two solitons, qf = q 2 , 
q^- = —q\ (and similarly for q^), so (6.7) simplifies to 

r7 = ^i^ ( /UlJ R)q 1 -q 1 , (6.9) 

as found in [4]. This is maximally attractive when qi = — qi, independent of the orientation 
of qi relative to n. Depending on the details of V, there may be higher order terms which 
break this symmetry. Note that the choice of eigenvectors e\,e 2 is purely arbitrary in this 
case (since /ij~ 2 J v is the identity map), so the notion of standard orientation is similarly a 
matter of convention. 

Given the above, it is interesting to determine what properties of V, v will enforce Hi — /jl 2 - 
Let G C 0(3) denote the subgroup of the isometry group of S 2 leaving V invariant, and 
G v C G the isotropy subgroup of v in G. There is an induced isometric action of G v on 
T V S 2 which commutes with J v , so the eigenspaces of J v are invariant under G v . Hence, if 
Hi 7^ /i 2 , T V S 2 must have a line (in fact, an orthogonal pair of lines) invariant under G v . 
Now G v C 0(2) (the isometry group of T V S 2 ), and an element of 0(2) fixes a line if and 
only if it has order 2. Hence, unless G v = 1 or Z 2 or Z 2 x Z 2 , T V S 2 has no such fixed line, 
and we conclude that fii = f^2- It does not immediately follow that U simplifies to (6.9) 
however, since this also requires qi = q 2 . The last condition follows if we assume (as is 
certainly plausible) that the single soliton is G v equivariant, as we will now demonstrate. 

So assume G v is nontrivial and different from 7L 2 and Z 2 xZ 2 . Having chosen eigenvectors 
€1,62 we have an induced isomorphism between the vector spaces R 2 (physical space) and 
T V S 2 , which we can use to transfer the G v action to E 2 . Then a map 4> : R 2 — > S 2 is G v 
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equivariant if g<f>(x) = 4>(gx) for all g G G v . Let X : R 2 — > T V S 2 be the asymptotic field 
defined in (6.2) so that (f> — v + X + • • •, namely (since Hi = /x 2 ), 

= d^itl!l {qixei + g2ye2 ). (6 .io) 
In H\r 

Equivariance of <fi implies <?X(x) = X(gx), and hence every g e G v commutes with 
diag(gi,g 2 ). But, if gi 7^ g 2 , this implies G v has only elements of order 2, which, by as- 
sumption, is false. Hence if the single soliton is G v equivariant, q\ = q 2 - 

The potential studied in the current paper, (2.4) with N > 3, nicely illustrates this 
symmetry analysis. If we choose boundary value v = (0, 0, 1) then G v = and /ii = /i 2 = 
to. Further, (at least for iV = 3) the single soliton is observed to be equivariant, so 
Qi = <?2 and the long range forces are as described by (6.9). By contrast, if we choose any 
of the other vacua, for example v = (1,0,0), then G v = Z 2 and there is no reason why ^ 
and /z 2 should be equal. The eigenvectors must be invariant under the Z 2 action, which fixes 
them, up to orientation, as (0, 1,0) and (0,0, 1). The corresponding eigenvalues are 2N 2 m 2 
and 0, so with our choice of conventions, /ii = 0, /i 2 = \/2Nm, €1 = (0, 0, 1), e 2 = (0, —1, 0), 
and the long range forces are as described in [11]. 

It would be interesting to see how this analysis generalizes to the (3 + l)-dimensional 
Skyrme model. 



Appendix B: the stress-energy of a baby Skyrme lattice 

In sections 4 and 5, baby Skyrme models with doubly periodic boundary conditions were 
studied. Equivalently, the model was put on a torus C/A, where A is a period lattice, chosen 
in this case to be rhombic, A* = {(n + me l7T / 3 )L : n,m E Z}, where the side length L > 
is determined numerically. Of course, given any lattice A, and any topological charge B, 
one would expect the baby Skyrme model on C/A to have an energy minimizer of charge B. 
Not all such doubly periodic solutions can be meaningfully interpreted as soliton lattices, 
however. For fixed B, we have a map which sends the lattice A to the energy of its charge 
B minimizer, and to be a genuine soliton lattice, A should be (at least a local) minimizer of 
this map. That is, the energy of a soliton lattice should be stationary under variations not 
just of the field 0, but also of the lattice A. This condition can be usefully reformulated in 
terms of the stress-energy of the field 0, as we now show. 

All tori C/A are diffeomorphic via real linear maps, so we can fix a standard lattice, 
for example A*, and consider every other torus C/A to be identified with C/A*, but with 
a nonstandard metric g (the pullback of the usual metric on C/A by the diffeomorphism 
C/A* — > C/A). So now the domain of 0, call it M, is fixed as a smooth manifold, but has 
a Riemannian metric g which varies as we vary A. In order to be a soliton lattice, a field 
: (M, g) — > S 2 should be a critical point of E under all smooth variations of 0, and all 
variations of g arising from changing A. This leads us to compute the variation of E(<f>,g) 
with respect to g. 
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It costs no effort to put the computation in a general geometric setting. So let M be 
a compact oriented n-manifold, (JV, h) be a Riemannian manifold, u be an n-form on N 
and V : N — > R be a smooth function. For a given metric g on M, the energy of a map 
: M — » iV is, by definition, 

W^) = ^l|d0|| 2 + i||0^|| 2 + / M ^(0)vol 5 (6.11) 

where || • || denotes £ 2 norm and vol s is the volume form on M associated with g. We wish 
to compute the variation of E(<f>,g) with respect to g. Let g t be a smooth curve in the 
space of Riemannian metrics, with go = g and dtgt\t=o — £■ Note that e is (like g) a section 
of 2 T*M, a real vector bundle over M of rank \n{n + 1). This bundle inherits a fibre 
metric from g, which we denote (•, •), defined as follows: let e±, . . . , e n be a local orthonormal 
coframe on (M, g); then we demand that <8> + ® e,) : « < j} is a local orthonormal 
frame for (© 2 T*M, (•, •)). The key fact is: 

Proposition Let gt be a smooth, one-parameter family of metrics on M with go = g and 
0t0t|t=o = £■ Then, for fixed : M ^ iV, 



= ^(5(0),^ 
t=0 Z 



where 



5(0) = (i|d0| 2 - i|0* w | 2 + V(0))s - 0*/i 



which, following the terminology of harmonic map theory, we call the stress-energy tensor. 
Note that, like 5 and e, S(<f>) is a section of 2 T*M. 

Proof: We compute separately the variations of the three terms in £, which we denote £"2, £"4 
and E respectively. The first term £ 2 = |||d0|| 2 is the Dirichlet energy of 0, whose variation 
with respect to g is [17] 



,=0 22 



W|d0| 2 £-0*M) i2 . (6.12) 



In the course of the proof of this, one finds that 



<9 t | t=0 vol 3t = ^tre vol p = ^(e,g)vol g . (6.13) 



It follows that the third term, E = J M V((j))vol g , has variation 



d 

-^E ((f),g t ) 



= l(V(<f>)g,e) L2 . (6.14) 
t=o z 



Finally, to handle the middle term £4, we must compute the variation of the Hodge iso- 
morphism * t : Vt n {M) — > Q°(M) defined by g t . Let /1 be a fixed n-form on M. Then, by 
definition * t fi = f t where \i = f t vol gt . Differentiating this with respect to t at t — yields 

= (d t f t )\ t=0 vo\ g + ^f (e,g)vo\ g (6.15) 
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using (6.13). Hence 



and so 

d 
dt 



t=0 



1 d 
2dt 



dt(*tn)\t=o = ~2^ £ ^) *oH 



[ <p*oj A * t (f)*uj = — - [ (e,g)4>*u A 
Jm 4 Jm 



(6.16) 



— a {e,\4>*u\ 2 9)l> (6.17) 



which completes the proof. 



x 2 / 



x, 



Figure 7: The period parallelogram of a generic torus. 



We now return to the problem of finding criticality constraints on a baby Skyrmion 
lattice. So from now on, M = T 2 = C/A, where A is some fixed lattice, = S 2 , oo = 
k x (area form on S 2 ) and V is some potential. Without loss of generality, we can assume 
that one of the periods of A is positive-real, and the other has argument a G (0,7r). Then 
the period parallelogram is as depicted in Figure 7. There is a three-parameter family 
of variations of A, up to isometries, generated by (a) homothety (uniformly scaling the 
parallelogram), (b) stretching the parallelogram horizontally, and (c) varying the interior 
angle of the parallelogram through a. We claim that the corresponding variations of the 
metric g = dx\ + dx\ are tangent to the symmetric bilinear forms 

£(a) = 9 
E( b) = dx\ 

£( c ) = — 2dx\d%2 + cot adx\. (6.18) 

Of these, £( a ) and £(&) are clear. To obtain £( c ), we note that varying a as a + 1 is equivalent 
to defining the inner product between the fixed pair of unit vectors X±, X2 to be cos(a + 1), 
while keeping their lengths unchanged (see Figure 7). Hence, 

e(X 1 ,X 1 ) = e(X 2 ,X 2 ) = 0, e(X 1 ,X 2 ) = e(X 2 ,X 1 ) = -sin a. (6.19) 

Now 

d d 1 

X 1 and t^ = - (X 2 -cosaXi), (6.20) 



dxi dx 2 sin a 
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so 



e{d/dx\, d/dxi) = 0, e(d/dx\, d/dx 2 ) = — 1, e(8/8x 2 , 8/8x2) = cot a 



(6.21) 



which is equivalent to (6.18). If <fi : M — > S 2 is a soliton lattice, then, by the Proposition, 
S((f>) must be L 2 orthogonal to each of £( a ), £(6)j £(c)- Hence, <S , (^>) must be L 2 orthogonal to 
any section of 2 T*M in the span of these, for example 



£( ) = 5, £( 6 ') = dx 1 dx 2 , 
Now ((?, g) = n = 2, and (0*/i, (?) = |d0| 2 , so 



£rr' 



(c') 



(5(0),£ (a) ) L2 = & j M {-\Wu\ 2 + V(<P)) vol, = 0, 

which coincides with identity (4.2). Further, (e^),g) = (e^),g) = 0, so 

(S{6),£(v\) — <^> {4>*h,dx 1 dx 2 ) — \ ^ — • — dx\ dx 2 — 
which is identity (4.4), and 



(6.22) 



(6.23) 



(6.24) 





8(f) 


2 


d<j> 






8x\ 




8x 2 


1 



dxi = (6.25) 



which is identity (4.3). 

We conclude by making two remarks. First, it is interesting that we get the same integral 
constraints on a baby Skyrmion lattice for all tori C/A. Second, in the case where the 
minimizer is holomorphic (e.g. the alternative D 3 model in section 5, or Ward's model 
[12]) only the scaling constraint (6.23) is nontrivial, since is conformal and 8/8xi, 8/8x 2 
are orthonormal, so 8(f) jdx\ and 8(p/8x 2 have equal length and are orthogonal, pointwise, 
and hence (6.24) and (6.25) hold automatically. 



Note added 

The preprint [18] contains several ideas related to those appearing in this paper. In particular, 
there is a similar proposal to identify quarks inside Skyrmions, and a suggestion that many 
new Skyrmions might be found as pieces of the Skyrme crystal; in the same way that the 
polyiamond solitons may be viewed as pieces of the soliton lattice. 
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